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ABSTRACT 

We have developed a new technique called Direct Shear Mapping (DSM) to measure 
gravitational leasing shear directly from observations of a single background source. 
The technique assumes the velocity map of an un-lensed, stably-rotating galaxy will be 
rotationally symmetric. Tensing distorts the velocity map making it asymmetric. The 
degree of leasing can be inferred by determining the transformation required to restore 
axisymmetry. This technique is in contrast to traditional weak leasing methods, which 
require averaging an ensemble of background galaxy ellipticity measurements, to ob¬ 
tain a single shear measurement. We have tested the efficacy of our fitting algorithm 
with a suite of systematic tests on simulated data. We demonstrate that we are in 
principle able to measure shears as small as 0.01. In practice, we have fitted for the 
shear in very low redshift (and hence un-lensed) velocity maps, and have obtained null 
result with an error of ±0.01. This high sensitivity results from analysing spatially re¬ 
solved spectroscopic images (i.e. 3D data cubes), including not just shape information 
(as in traditional weak leasing measurements) but velocity information as well. Spirals 
and rotating ellipticals are ideal targets for this new technique. Data from any large 
IFU or radio telescope is suitable, or indeed any instrument with spatially resolved 
spectroscopy such as SAMI, ALMA, HETDEX and SKA. 
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1 INTRODUCTION 


Weak gravitational leasing maps matter distributions in the 
universe, both baryonic and dark (e.g. [Kaiser &: Squires] 
119931). This paper explores the enhanced potential of weak 
leasing in three dimensions ( Morales|2006 Blain|2d02 1. We 
describe a methodology to obtain a shear measurement from 
a 3D data cube of a single weakly lensed galaxy. This tech¬ 
nique will allow us to measure the size and shape of dark 
matter distributions around individual galaxies at low red- 
shifts. 

Conventional 2D weak leasing techniques rely on measur¬ 
ing the (two-dimensional) shapes of many (> 100) images 
in a field. With the assumption that the images in the field 
should have no preferred orientation, one can infer the pres¬ 
ence of a shear field if any correlation in alignments is de¬ 


tected. There are a number of statistical approaches (Bridle 
et al.||2010| [Reymans et al. 112006} [Kaiser et al.[[T995[ [Re 


fregier & Bacon 20031 used to perform this analysis. The 


statistical uncertainty of a particular weak lensing survey 
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depends on the total survey area and the number density 
of perfectly measured galaxies. The dominant source of un¬ 
certainty in all of these methods is shape noise: the error in 
the measurement due to the intrinsic and random orienta¬ 
tions of the images in the sample. Using these techniques, 
one would typically need ~ 10 objects to measure a shear of 
10% on arcminute scales. 


A powerful enhancement of the weak lensing method was 
proposed by Blain (20021, followed by Morales (20061. The¬ 
oretically, the rotation curves of regularly rotating elliptical 
and spiral galaxies will have maximum and minimum val¬ 
ues in the projected velocity maps. These coincide with the 
major and minor axes of the projected 2D image. In gravita¬ 
tional lensing, the Equivalence Principle requires that pho¬ 
tons of different energies are affected similarly. Thus weak 
lensing will shear the velocity maps and distort the 2D im¬ 
age, causing the angle between maximum and minimum ro¬ 
tation axes to deviate from 90°. While there are other kine¬ 
matic effects which also introduce perturbations into the 
velocity map of a galaxy, weak gravitational lensing has a 
unique signature. An un-sheared velocity map is symmet¬ 
rical about the major and minor axis, while a sheared ve¬ 
locity map loses these symmetries. Since we can predict the 
properties expected in an un-sheared velocity map (namely 
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that the angle between the major and minor axis is orthog¬ 
onal, or put another way, that the velocity is symmetrical 
about these axes), a metric which measures the presence 
and strength of the shear field between the observer and the 
background galaxy can be designed. 


Blain (20021 suggested using the distortion in the rotation 


curve to measure shear. He fitted for shear in an inclined 
ring model using a Monte-Carlo routine, concluding that 
while currently a typical shear could not be measured by 
this method, it would be possible to measure it with fu¬ 
ture higher resolution surveys. Morales (20061 suggested 
a similar method, but fitted for the entire velocity map. 
This has the advantage of utilising additional information 
from the velocity map, and potentially avoiding problems 
the concentric-ring method would encounter when htting 
for shear in warped or disturbed disks (since in that case 
each concentric ring will recover a different shear). Morales’ 
method involved measuring the angle between the major 
and minor rotation axes, and obtaining a measure of the 
shear strength from the deviation of the fitted angle from 
orthogonal. 

We have extended Morales’ work, and developed a technique 
to measure the shear vector directly from the 3D data cube 
of a single weakly lensed galaxy. Our technique is called 
Direct Shear Mapping (DSM) and utilises a Monte-Carlo 
Markov Chain (MCMC) algorithm to search for asymme¬ 
tries in the 3D data. The MCMC fu nction used is called 
emceeQ | |Foreman-Mackey et al.|[2013 1. We have performed 
a suite of tests on simulated data to characterise the effi¬ 
ciency and accuracy of the fitting algorithm, and to under¬ 
stand the limits of its fitting range. Additionally we have 
tested the technique on un-sheared data to ensure we can 
recover a null result at low redshift, and determine realistic 
systematic errors. 

Using a symmetry search method rather than focusing on 
the velocity axes directly allows us to use all the available 
data, giving a statistically better fit. This method is unique 
among weak lensing measurement methods in that it uses 
the velocity map of an object to measure the shear field, 
and does not rely on fitting for the shape of the galaxy. 
Rather than obtaining a single measure of the shear over a 
wide area of sky, we can determine a shear value for a single 
background galaxy. Hence we are able to measure the mass 
of an individual foreground dark matter halo. 

We have demonstrated that we can measure shears as 
small as 7 ~ 0.01 in realistic simulations and expect to 
extend this to observational datasets (Taylor et al, in 
prep.). Two-dimensional weak lensing methods are able to 
statistically measure the mass and structure of an individual 
dark matter halo and establish its relationship to other 
observables, for example baryonic mass (eg|Velander et alT 


120141). It may also be possible to probe for significant 
substructure in the shear field on small angular scales 
by making a number of independent shear measurements 
around a dark matter halo, a problem which shape-htting 
and moment-measurement methods cannot address, since 
they are blind to variations on smaller scales. We already 
have a first probable detection (Taylor et al, in prep.), and 
there are a number of existing and upcoming instruments 


^ the Emcee Hammer, http://dan.iel.fm/emcee/ 


and surveys which will produce data products well suited 
to our technique, including observations with SPIRAL 
I Kenworthy et al. 2001| , SAMI (Groom et al. 20121 and 
CALIFA ( Sanchez et al.|[2012 l in the optical, and ASKAP 


and the SKA at radio frequencies (Johnston et al. 2008 
Schneider|2000 Braun|2007). 


This paper is organised as follows. In section we briefly 
review weak lensing theory. In section we describe DSM, 
and the structure of the DSM fitting algorithm. In section]^ 
we describe the generation of synthetic galaxy data, and 
present the results of a suite of systematic tests of the DSM 
method. In sectionj^we present fits to low-redshift data from 
the literature, demonstrating the ability of DSM to recover 
a null result. Conclusions are presented in section 


2 WEAK LENSING 


Lensing theory is well developed, and there are many deriva- 


tions of the relevant equations (e.g. 

Mortlock|| 1999 Bartel- 

|mann & Schneider|2001| |Bartelman 

20Tof^ 


Light travels along geodesics, and in the presence of massive 
bodies, geodesics are curved. As light passes massive bod¬ 
ies, its path is deflected, and a background image will be 
distorted, magnified, shifted and duplicated. Weak gravita¬ 
tional lensing is the term used to describe minimal distortion 
of the background source with no observable multiplication 
of the source image. 

For this paper the focus is on linearisable weak lensing, 
that is, lensing where the distortions are small, and do 
not vary across the source. This class of lensing can be 
expressed as a single transformation matrix. It is this 
feature which forms the basis of the method by which the 
DSM shear values are determined. We will also restrict 
ourselves to individual lenses (e.g. galaxies or clusters), 
however this method is equally valid for measuring the 
cosmic lensing signal. 

If one assumes the length scales of the lensing mass distribu¬ 
tion are much smaller than the observer-lens distance, then 
the thin lens approximation can be used. In the thin lens ap¬ 
proximation we project the lensing mass distribution onto a 
plane perpendicular to the observer’s line of sight to obtain 
the surface mass density: 

m) = jdzp{iz), ( 1 ) 

where ^ is a set of coordinates in the lens plane, 

and z is the third coordinate parallel to the observer’s line 
of sight. Using the weak field approximation, the deflection 
angle of light passing the lensing mass is then given by 



where E(^) is the (projected) surface mass density and 

M{< I) = I d"^'E(f) (3) 

is the mass enclosed within The surface mass density can 
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be written in a dimensionless scaled form k(^) = E(^)/Scr-, 
where 

-(4) 

4nGDdDds ^ ’ 

is the critical surface mass density. E ^ Ecr yields multiple 
images (i.e. strong lensing) while E ^ Ecr gives only one 
(i.e. flexion and weak lensing). 


The lens equation is 

/3 = e (5) 

where 8 = ^/Dd is the apparent angular separation of the 
lens and source, (3 = 8/Da is the true angular separation 
of the lens and source and r] is the projected distance be¬ 
tween the lens and source in the source plane. Dds is the 
angular diameter distance between the lens and source, 
is the angular diameter distance between the observer and 
source and Dd is the angular diameter distance between the 
observer and lens. 

In the thin lens limit, if the bend angle is slowly varying with 
position, then the lens mapping can be locally linearized. 
Equation Q can then be re-expressed as a coordinate map¬ 
ping between the lensed and un-lensed coordinate systems 

P = Ae, ( 6 ) 


where A is the Jacobian of transformation and is given by 


A{e) 


dS ^ dOidOj 

/ 1 — K — 7l —72 
\ -72 1 - K -I- 71 


( 7 ) 


where j3 is the angle to the true source location, k is the 
convergence (as defined previously), and 71 and 72 are the 
shear components given by 


2^881 88^' 88^882' 


( 8 ) 


and 7 = ( 71 , 72 ). The magnification is given by /r = = 

(i_B,)2_[.y|2 ■ A is the most general form of the linearised 
weak lensing transformation, but not the most useful for our 
purposes. In the weak regime, k has only a small effect as 
E <C Ecr-. Furthermore, since the DSM algorithm compares 
the velocity map of a galaxy to reflections of itself, and k 
produces a spatial magnification only, it will have no effect 
on the fit. With this in mind, the k term in Equation 0 
is discarded for simplicity, so that we can define a reduced 
transformation 




1 — 71 —72 \ 

-72 1 + 7i / ’ 


( 9 ) 


and 71 = I 7 I cos( 20 ) and 72 = | 7 |sin( 20 ). Ar can be ex¬ 
pressed in a simpler form by selecting a coordinate system 
such that the shear vector is at a fixed angle with respect 
to the coordinate axes ( Morales|2006 |. We choose to define 
our coordinates such that the shear vector lies at 45° to the 
j/-axis, to obtain 


7 I 45 = 


1 

-7 



( 10 ) 
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where 7 = | 7 |. It is an important property of 7 I 45 that it is 
invertible and unitary, i.e. 


Aa^ — 


( 11 ) 


1 — 7 ^ 

and 71 ^ 5^7145 = ^ 457145 ^ = I. To recover the full shear vector 
in a generalised coordinate system with angle 8t to the shear 
vector, from a measured shear of 7meas = I 7 I, one simply 
projects onto the new coordinate system: 


7 = ( 7 I 172 ) = (7measCOs(26't),7measSin(26lt)). (12) 


The angle 8t can be obtained by knowing the direction to 
the lensing mass. In the cases considered in this paper the 
lensed galaxy is at low redshift and so the direction to the 
lensing mass can be determined. We now have a linearised 
lensing transformation matrix 7 I 45 with which to map our 
image coordinates according to a shear strength 7 . 


To set up the problem, we first define terminology and the 
relevant coordinates. There are four relevant coordinate sys¬ 
tems to consider: 


(i) The intrinsic galaxy coordinates in the source plane, 
prior to lensing: (m, n). 

(ii) The observed galaxy coordinates in the detector 
plane: (m',n'). The prime notation refers to the observed 
frame, in which a shear signal is present. 

(iii) The detector frame coordinates: {x',y'). These coor¬ 
dinates correspond to the pixel coordinates of the detector. 
Note that since they sit in the plane in which a shear signal 
is present, they are also primed. 

(iv) The intrinsic frame resulting from taking the detector 
coordinates and projecting them back into the source plane: 
ix,y)- 

While it may at first seem counter-intuitive to use primed 
notation to denote the observed frame, it is natural to do so 
considering it is the observed frame that it is lensed, and the 
intrinsic source plane which is not lensed. Note that at no 
point are the coordinates of the lensing galaxy referred to. 
The only coordinate systems used are those of the intrinsic 
plane of the source galaxy, and the plane of the detector. 
The detector coordinates are centred on the middle of each 
pixel, and relate to the un-lensed source plane coordinates 
(x, y) through 


= Ai5 


(13) 


The detector plane galaxy coordinates are the coordinates 
of the galaxy as it is observed on the sky, and are related to 
the intrinsic source plane galaxy coordinates through 


m! 

/ 

= 7145 

m 

n 


n 


(14) 


where the ‘source plane’ refers to the plane in which the 
true galaxy sits, (m, n) will be aligned at some angle cp with 
respect to the source plane, and will have a different origin 
to the {x, y) coordinates system, i.e. 


X — Xo 

_ y-yo 

where TZ((p) is a rotation matrix, and 


Xz 


X — Xo 



y-yo 


= 7^(<^) 


m 

n 


(15) 


(16) 
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where (xo,yo) are some translation of coordinates, such that 
in the (xz,yz) coordinate system the galaxy is centred on the 
origin. This relationship between {x„,yT) and (m, n) will be¬ 
come relevant later when discussing the fitting methods em¬ 
ployed by DSM. 

Finally, since .445 is unitary, to undo a leasing transforma¬ 
tion, one simply needs to apply the inverse matrix to 
the lensed coordinates; 


-^45 


X„ 

y'z 


— -^45 -445 


Xz 


N 

H 



. y^ 


(17) 


3 METHOD 

3.1 The DSM symmetry method: 


The aim of the DSM algorithm is to use an MCMC htting 
algorithm to find the best fit shear value in an observed 
object, taking as its input a velocity map (V), an error map 
(E) and a mask (M). This is done by looking for axes of 
symmetry in the velocity map. As shown in Equation (121, 
when combined with knowledge of the direction to the lens, 
the best fit shear value provides a measurement of the total 
shear field along the line-of-sight to the observed source. A 
set of such measurements can then be used in a number of 
scientific applications. The lensing object can be either a 
cluster or field galaxy. 


The key idea with this technique is that the projected 
velocity map of a rotationally supported galaxy is sym¬ 
metric about two axes (the axis of maximum and of zero 
projected rotational velocity). We are specihcally looking at 
low redshift background sources, so that we have sufficient 
angular resolution to apply this algorithm. The low redshift 
of the source and lens means the direction to the lensing 
mass is well determined, assuming the foreground luminous 
galaxy is at the centre of mass of its dark matter halo. 
Thus we measure the mass of the dark matter halo of the 
foreground galaxy, through weak shear and possibly flexion. 


For circular orbits (such as in spiral galaxies), the projected 
velocity onto the plane of the observer is (u(r)m/r) sin(^), 
where r is the distance from the centre of the galaxy, m is the 
projected distance along the major axis and ^ is the angle 
between the line-of-sight and the galaxy’s rotation axis. This 
results in the maximum rotational velocity along the m axis, 
and the minimum along the n axis, with the angle between 
them being 90°, where m and n are the projected distances 
along the major and minor axes of the galaxy respectively. 
A diagram of this geometry is presented in Morales (2006 \, 
Figure 1. 

This orthogonality implies that we should find three symme¬ 
tries in the velocity map of a rotating galaxy: symmetrical 
about the major axis, anti-symmetrical about the minor axis 
and anti-symmetrical about the diagonal between the two. 
Once a galaxy is lensed however, the shape of the galaxy 
is distorted without the frequency of the light received 
at each pixel being changed. Since the rotational velocity 
of a galaxy is determined via the red- or blue-shifting of 
a known emission frequency, it follows that the apparent 
rotational velocity of the galaxy will not change under 
lensing, but the shape will. In other words the image will 


be distorted along an axis in the direction of the shear 
field, so that the velocity axes m, n will no longer have an 
angle of 90° between them, resulting in the loss of all three 
symmetries from the image. DSM uses this deviation from 
symmetry to fit for the shear field that the galaxy light has 
passed through. For further reading on the topic of velocity 
maps in rotationally supported galaxies see |van der Kruit| 


& Allen (19781. 


DSM fits for the gravitational lensing shear field which best 
restores symmetry to the image. Since in the weak lens¬ 
ing regime the light distortion is small, it can be linearised 
and expressed as a single transformation matrix. Thus if a 
transformation matrix can be found that restores symmetry 
to the image, this can be inverted to End the transformation 
that resulted in the initial distortion, and hence the shear 
strength. 

To measure the degree of asymmetry in the image, the DSM 
algorithm uses an MCMC routine with many walkers, which 
aims to maximise a likelihood function by stepping through 
the parameter space of possible solutions for (x' ,y'), R and 
S (the definitions of which are given below). To do this, the 
MCMC routine iterates over the following steps: 

(i) Translate the image plane galaxy coordinates (m', n'), 
so that the galaxy centre sits at the origin in the image plane 
coordinate system (x',y'). 

(ii) Apply the following coordinate transformations to a 
regular coordinate grid, representing the coordinates of the 
pixels in a (possibly sheared) velocity map image: 

S"^ =:> => T ^ K =» S, (18) 

where S (S~^) is a shear (inverse-shear), R (R“^) is a rota¬ 
tion (inverse-rotation), and T is a reflection about one of the 
three symmetries in the image. The coordinate set resulting 
from these transformations is used to map the input image 
to three reflected trial images. 

(iii) In the detector coordinates (a;', y'), compute the per- 
pixel differences between the reflected and reference (i.e. 
original) images, resulting in three ‘difference images’. 

(iv) For each difference image, compute the per-pixel 
(giving three sets of values) 

(v) For each set of per-pixel x^’s, compute the likelihood 
value. 

(vi) Add the three sets of per-pixel likelihoods, and then 
sum over the pixels to obtain the total likelihood for this fit. 
The shear is obtained from the MCMC maximum likelihood 
fit. 

A diagram of steps (i) to (iii) is shown in Figure 


The structure of the DSM algorithm is now described in 
detail. 


3.2 The fitting algorithm 

The core of the fitting algorithm, trial_fit, is struc¬ 
tured such that each set of trial parameters is input into 
TRIAL_FIT, and a log-likelihood is returned. The set of 
parameters that results in the maximum likelihood is found 
by inputting trial_fit into a Monte-Carlo routine. The 
module emcee is used in preference over other modules 
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Transformations of the observed image projected in 2 spatial dimensions 




(h) 



n = Detector plane (blue) X 
□ = Source plane (red) 




(f) 


Y 




Figure 1. A schematic of the transformations in the fitting process. Blue (solid lines) corresponds to the object in the detector plane. 
Red (dashed lines) indicates the object projected back into the source plane. It is important to understand that while each step is shown 
individually here, this is for illustrative purposes, and in the DSM algorithm all of the steps shown are done in one transformation. This 
is in order to minimise interpolation errors. Since the object is initially in the detector plane, and ends up in the same plane, so long 
as all of the transformations shown here are done in one step (in the order shown), the only point at which interpolation is necessary 
is in the final step, when computing the residual. The steps are as follows: (a) is a translation to the origin of the detector coordinate 
system {x\y') from the detector plane galaxy coordinates (m',n'). (b) is an inverse-shear, taking the object from the detector coordinate 
system to the un-lensed source plane equivalent, [x^y). (c) is a rotation, the aim of which is to leave the object aligned with the axes 
of the detector coordinates (x',y'). In (d) reflections are made about the x and y axes of the [x'^y') coordinate system, (e) reverses the 
rotation applied in (c). (f) re-applies the shear removed in (b), taking the galaxy back to the detector plane, (g) is a translation, taking 
the galaxy back to its original location, so that it is centred on the (m',n') coordinate system, (h) is the computation of the residual, 
and interpolation back onto the detector coordinates [x'^y'). The combination of transformations in this series is clearly not the correct 
set of values, as there is a nonzero residual apparent in (h). 
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because it employs an ensemble of walkers. This is advan¬ 
tageous for this application since the likelihood surface 
features local minima (representing incorrect solutions, 
expected to be a result of discretisation and interpolation 
errors). The Monte-Carlo routine is invoked by a function 
RUN_EMCEE. The order of processes in run_emcee is as 
follows: 

There are five parameters the algorithms seeks to optimally 
measure; shear ( 7 ), galaxy position angle {9), galaxy centre 
in the {x, y) coordinate system (xo and j/o), and galaxy cen¬ 
tral velocity (wo). A set of initial guesses for the parameters 
are specified, followed by the parameters for emcee to run 
(i.e. number of walkers, number of initial steps to discard 
(‘burn-in steps’), and total number of iterations). A set of 
initial walker locations is then generated in a ‘ball’ around 
the initial guess values. 

A log probability function is defined, which is passed to the 
htting parameters. The log probability function hrst deter¬ 
mines whether all the parameters to be htted lie within a 
pre-set range of realistic possible values. 

If they do, it calls trial_fit. trial_fit returns the total 
likelihood and the total ■ The burn-in steps are run, and 
the sampler is then reset so that the burn-in steps aren’t 
counted toward the final ht. emcee is then run for the full 
number of iterations. run_emcee then returns the object 
‘sampler’ and the array of acceptable walkers. The fitted 
values of the parameters are obtained by taking the mean 
of the fitted values of the walkers, then taking the mean of 
the walkers (in each dimension, for each parameter) 


XQ = {{X0i)j), 

Vo = {{yoi)j), 

Vo = {{voi)j), (19) 

7 = ((7i)j), 

ft = mA 

As mentioned above, the core of the algorithm is trial_fit, 
which generates the likelihoods, and so some time is taken 
here to discuss its form. trial_fit takes as its input the set 
of parameters to test, and the data in the form of a velocity 
map {V), an error map {E) and a mask (M). 


Three reflection matrices are then defined, for reflection of 
the galaxy coordinates about a major axis, a minor axis 
perpendicular to the major axis, and about both axes. If the 
line of reflection is hxed to pass through the origin, and has 
an angle with respect to the x-axis of uj, then the generalised 
reflection matrix is 


f cos( 2 a;) sin( 2 a;) \ 
y sin( 2 w) — cos( 2 a;) ) 


( 20 ) 


DSM is structured such that a; = 0, so that the reflec¬ 
tion axes always align with the detector coordinates {xAvA- 
With oj = 0, the matrices to reflect the coordinates about 
the major axis, the minor axis and both axes are: 


• About the major axis: 



( 21 ) 


• About the minor axis: 


Ta 


• About both axes: 



T 3 = T 1 T 2 



( 22 ) 


(23) 


where the subscripts denote the different axes of reflection 
(1 = major, 2 = minor and 3 = both). 


The inverse shear and rotation matrices to be tested are 
then defined as 


S ^ = (1 - Xtrial) ^ 


1 

'Jtrial 


'Atrial \ 

1 )' 


(24) 




COsixf^trial ) 
Sin(<()fT.iai) 


COs{4>trial) J ’ 


(25) 


followed by the corresponding shear and rotation matrices 


S = 


1 

~'ytrial 


'Atrial 

1 


(26) 


f COs( (ptrial) 
y Sin( (j^triaA 

k COsi^cl^trial) 

y sin((()t7.jai) 


- Sm{-(j)trial) 
COs( (j)trial ) 

Siri{(j)trial) \ 
COs{4>trial) J ’ 


(27) 


where Xtriai is a trial shear strength and (p-triai is a trial 
value for the position angle of the object in the source plane. 


For each reflection (about the major axis, the minor axis and 
both), the matrices are multiplied together in the order: 

^1,2,3 = S ffi Ti 7,3 R (28) 


such that an inverse shear of Xtriai is the hrst operation to 
take place, followed by a rotation by —cptriai, and so on. This 
series of matrix transformations is then applied to a set of 
x'z and 2/2 coordinate arrays (i.e. the detector coordinates) 
to produce trial ht coordinates; 


^21,2,3 ^1,2,3 Xz, 

3 / 21 , 2,3 = Ti,2,3 y'z, 


(29) 


which are then translated back to the {x, y) coordinate sys¬ 
tem 


* 1 , 2,3 — * 21 , 2,3 + * 0 , 
j/ 1 , 2,3 = yt 1 , 2,3 + yo- 


(30) 


Since these coordinates have been transformed, they will not 
necessarily be integer values and so it is necessary to inter¬ 
polate the data when drawing values from the input images 
with the trial coordinates (* 1 , 2 , 353 / 1 , 2 , 3 )- In th® DSM algo¬ 
rithm bilinear interpolation is used. To understand the bi¬ 
linear interpolation process, it is useful to think of the trans¬ 
formed coordinates as a distorted grid, which we are trying 
to project back onto a regular grid. In this case each dis¬ 
torted grid square will fall across four regular grid squares, 
or alternatively each regular grid square will have part of 
four different distorted pixels fall into it. In order to compute 
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the value of each regular grid square, we need to interpolate 
between the four pixel values it contains. We can work out 
what fraction of each distorted grid square I{x^,y^) falls into 
our regular grid squares V{p,q), and assign a value to the 
regular grid square which is equal to 

V{p,q) = Dtfytl{x\y^) 

+ fx<:fyt + lID yV + 1 ) 

+ fx<: + lfy^y + lI{x +1,2/ +1), 


where /^.t fyt denotes the fraction of the distorted grid square 
I{x^,y^) with coordinate {x^,y^) that falls within the regu¬ 
lar grid square V{p, q). In order to minimise interpolation er¬ 
rors, it is desirable to interpolate as few times as possible. To 
this end, the trial data is kept in a fractional form until the 
last step in the fitting process, the computation of the per- 
pixel X values, so that it is only interpolated once. In other 
words I{x*,y*), I{x^ + l,y*), I{x^,y*^ + 1) and I{x* + l,y^ + 1) 
are kept separate and transformed separately. This means 
that from a single input set of 3 data arrays (velocity map, 
error map and mask), we have 3x3x4 arrays until they 
are interpolated (3 input images, three reflections and four 
fractions). For simplicity therefore, in this paper the four 
fractional arrays (prior to interpolation) are written as a 
single array, since they are all treated identically up to the 
point of interpolation. 

With this in mind, the transformed (fractional) coordinates 
are used to create transformed data arrays 

^1,2,3 = ^(®l,2,3i 2/1,2,3)) 

El, 2 ,3 = -£'( 3 / 1 , 2 , 352 / 1 , 2 , 3)1 


where V*' is the transformed velocity map, is the trans¬ 
formed error map, and M* is the transformed mask. 

A set of three ‘difference images’ are then computed by sub¬ 
tracting the reference velocity map (i.e. original un-lensed 
image) from each of the transformed maps, along with an 
additive shift which is a function of vo, the form of which is 
dependent on which axis of symmetry is being used: 


AVi,2,3 = 


Vt - V 

2Vo — V2,3 — V, 


(33) 


where the two forms are due to the fact that the sign of the 
velocity map switches across the y-axis, but stays the same 
across the x-axis (thus the sign is reversed for reflections 
about the minor axis and about both axes, but remains the 
same for reflections about the major axis). 


The error associated with the difference images is given by 

A-Bi, 2,3 = 7(£?‘)l2,3+^"- (34) 

At this point any pixels in the velocity or error maps with 
bad or missing data are flagged. Any pixels on the edge 
of the object, where pixels containing velocity data would 
be interpolated with non-finite values are also flagged. A 
per-pixel x value is then computed for each of the three 
reflections 


Xl,2,3 — (AVi,2,3/Ai?i,2,3), (35) 

It is at this point that it is most advantageous to interpolate 
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the data, since the next step is to square the data. Interpo¬ 
lating prior to squaring the data reduces the interpolation 
errors since they will be linear, not squared. The interpo¬ 
lated X values are clipped to amax, a specified maximum 
allowed value. Any flagged data points are also penalised 
and set to a max- The setting of a maximum allowed value 
(Xmax avoids the situation where a few bad pixels with a 
very large error dominate the fit. Since the shear signal is so 
small compared to the errors from inaccurate data around 
the edges of the velocity map, if the x values are not clipped 
the fit would be dominated by minimisation of the residuals 
in the edge pixels. A check is made by eye to determine the 
location of the bad (clipped) pixels. The three interpolated 
X arrays are then added together and squared to obtain the 
per-pixel x^ value 

X =Xi+X2+X3- (36) 


The three error arrays (Ai?i, 2 , 3 ) also need to be interpolated 
at this point, since they are now squared in order to compute 
the per-pixel log-likelihood. The likelihood of a given pixel 
having a given value, assuming gaussian errors, is 


L = V —== 



(37) 


so the masked per-pixel log-likelihood is 

log 1/ = — ^ [ X^ + 2 log(Ai?i Ai52 AFis) + 3 log(27r) ] x M. 

(38) 

Note that the extra error terms arise due to the normalisa¬ 
tion in Equation ( |37[ and the log-likelihood has been mul¬ 
tiplied by the input mask M to ensure the number of pix¬ 
els included in the final log-likelihood value is the same for 
every trial fit. The per-pixel log likelihood, logL, is then 
summed over the pixels to obtain the total log-likelihood, 
logiC. TRIAL_FIT then returns log£, which is used by the log 
probability function logp. Clearly the smallest residuals in 
AVi, 2,3 will come from trial values which leave the object 
symmetrical about the three axes of reflection, i.e. the trial 
values corresponding to the true position angle and shear in 
the input image. 


4 MODELLING AND TESTING OF METHOD 

To fully characterise the sensitivity of the method, and the 
efficacy with which it recovers the true shear in an image, a 
suite of tests were performed on synthetic data. The galaxies 
were modelled using the following steps: 

(i) Generate a model galaxy, with a lensing signal present; 

(ii) Degrade the velocity and spatial resolution of the im¬ 
age (this is done after the lensing signal is imposed since 
these properties are chiefly due to seeing and instrumental 
uncertainties); 

(iii) Perform an MCMC implementation of the 
symmetry-search process, including input trial inverse- 
shears and rotations, and computing the x^ ^nd likelihood 
of the trial fit; 

(iv) Repeat the above processes systematically for a range 
of input parameter values; 

(v) Assess the robustness of the fitting process as a func¬ 
tion of each of the input parameters. 
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4.1 Generating the model data 

The model-generating algorithm, make_data, can be 
roughly broken into two steps; generating the clean maps, 
and then adding seeing and velocity measurement error. In 
order to add the seeing, it is necessary to subsample the 
data. To achieve this, the data is initially generated to have 
M X N X n pixels, where n is the number of subsamples per 
larger pixel. 

The data reduction process for radio data results in a 
single error value for every pixel in the velocity error map 
( [Thompson et aT]|2001[ ). In contrast, optical IFU data is 
constructed from many individually-measured spectra and 
so has a different error for every pixel in the velocity map. 
Since the error maps generated for the model galaxies are 
constant across the map, they are more analogous to radio 
data than to optical data. More sophisticated error maps 
which vary across the galaxy would be required to model 
optical data. This is left for future publications. 


To begin with, empty arrays are created for the veloc¬ 
ity map, intensity map, and mask. Arrays of coordinates 
are then created, running from zero to (xp, yp), with size 
(xp X ^yn, Yp X V^))i where Xp and yp are the lengths of 
the X and y axes of the final images to be made, and n the 
number of subsamples per pixel to be used when convolving 
with a Gaussian. These arrays are called Xi and yi. These 
coordinates are then inclined about the x-axis 


J/inci = (j/i)sin(j 7 ), (39) 

where t] is the inclination angle. These inclined coordinate 
arrays are then lensed and rotated by the required values 

Xi — /C Xi ^ 

Vinci — ^ Vinci ? 


(40) 


where 


K. = 


1 -7 \ 

f cos 9 

— sin 9 \ 

V -7 1 J 

— sin 0 

cos6' J 


(41) 


These coordinates are then used to generate a galaxy im¬ 


age with an exponential surface brightness prohle (Freeman 
1970 Binney fc Tremaine||2008 1 


I{x,y) = ioexp{-2R), (42) 

which gives the following expression for the velocity map 

V{x,y) = voR^ [Io{R)Ko{R) - hiR)KRR)] 

x' (43) 

X cos[arctan( R )], 

Vinci 

where io and vo are scaling constants which are set to unity, 
Ko, Ki and 7o, I\ are Modified Bessel Functions of the hrst 
and second kind respectively, and 


R = 


\J + y'in 

2Rd 


(44) 


where Rd is the disk scale radius. 


This inclination, lensing and rotation procedure assumes the 
galaxy has no thickness, which is a reasonable approxima¬ 
tion since the characteristic scale length of a spiral galaxy 
disc is typically several orders of magnitude larger than 


the thickness or scale height ( Gaensler et al.|[2008 l. Future 
implementations of the algorithm will investigate this 
assumption, and if the effects of thickness are non-negligible 
it can be incorporated into the modelling. Thickness will be 
more significant when considering elliptical systems, and if 
appreciable, will smooth out the shear signal. 

A mask is then created with the radius of the mask being 
set at some multiple of the disk scale radius, Rd- This mask 
ensures that only the uniformly-rotating inner parts of the 
galaxy are used in the fitting algorithm, and the outer 
regions, which are low signal-to-noise and often disturbed, 
do not bias the measurement. For real galaxies, the mask 
size is determined by inspection, including the regular 
inner regions of the galaxy. It is noted that there is no 
real kinematic signature that can mimic the lensing shear 
signal. 


The effects of seeing and imperfect spatial resolution on 
the velocity and intensity maps were simulated by convolv¬ 
ing the flux-weighted velocity map with a Gaussian kernel 
of width s which when combined with assumptions about 
model redshift and size, can be expressed in pixels. The 
Gaussian kernel is given by 


G{x,y,s) = 


2'ks^ 


exp 


-( 


2 I 2 

X +y , 
2s2 • 


(45) 


so that the flux-weighted degraded image is given by 


V-w{x,y) = [V{x,y) x I{x,y)] (g) G{x,y,s). (46) 

The degraded images are then re-binned to the final size (i.e. 
the array sizes are re-binned from (xp x y^ x ^/n) to (xp, 
yp)), and normalised for the flux-weighting. 


Vdcg{x,y) = Vn,{x,y)/I{x,y) (47) 


Finally, the images have a measurement error (i.e. noise) 
added to the pixel values. The noise is drawn from a Gaus¬ 
sian distribution peaked around the central velocity, with a 
width specified as a fraction of the maximum velocity range 


Vdeg{,X, y) — Vdeg{x, y) -(- A^(Vo, V/rac), 
Vf rac — SVIVmax, 


where 


N{Vo,Vfrac) 


1 

rac 


exp 


1, Vo 

2 Vfrac 



(48) 


(49) 


where Vdcg is the 3D image with error in the velocity added, 
5V is the error in the velocity, Vfrac is the fractional error 
in the velocity, Vq is the central velocity, and Vmax is the 
maximum velocity range in the model. 


4.2 Sensitivity tests 

In order to characterise the sensitivity of the DSM algo¬ 
rithm, an ensemble of fits was computed for an exhaustive 
range of values of the input variables, specified below. The 
systematic tests were performed on synthetic data generated 
using MAKE_DATA, with many realisations of the data fitted, 
so that while each model had the same measurement er¬ 
ror characteristics, the error in each realisation was random. 
The sensitivity to the following parameters was investigated: 
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Figure 2. The error in the recovered shear values for an esemble 
of model realisations, as a function of the true input shear. The 
parameters used in the galaxy models are p = 7.5 pixels, di^cl = 
35°, Vmaxl&V = 100. The blue boxes with error bars show the 
mean and standard deviation in each bin, while the red circles 
show each individual realisation. There are 40 realisations per 
bin, with 12 bins over a range 7 = (—0.3,0.3). 


(i) Velocity signal-to-noise {Vmax/SV) in the range 
Vmax/5V = (2,1000). This incorporates the effects on the 
accuracy of the measured velocity from both spectral reso¬ 
lution and galaxy inclination angle; 

(ii) Sampling resolution (i.e. number of pixels across the 
galaxy) in the range p = (2,15). This is a proxy for galaxy 
angular size or redshift; 

(iii) Shear strength in the range of 7 = (—0.3, 0.3), where 
this is the same 7 as defined in Section to demonstrate 
the stability of the algorithm for strong shears; 

(iv) Inclination angle <j) in the range (j) = (10°,80°). This 
investigates the effect of a reduced number of pixels and 
a change in the projected shape of the velocity map with 
inclination angle. It does not investigate the effect of incli¬ 
nation angle on the accuracy of the velocity; this is taken 
into account in the investigation of velocity signal-to-noise, 
discussed above. 

(v) Once the other parameters had been investigated, see¬ 
ing (s) in the range s = ( 1 , 10 ); 

Initially, a set of 40 realisations were fit for a set of ‘fiducial’ 
values of the input parameters, to establish a baseline sen¬ 
sitivity. The values chosen were for a 20 kpc galaxy, masked 
beyond 1.5 Re at a redshift of 2 = 0.1, which gives a mask 
with a diameter of 15 pixels, and a galaxy with an Re of 5 
pixels. The fiducial inclination angle is Oinci = 35°, with a 
velocity signal-to-noise of VmaxISV = 100 . 

Once the baseline was established, the behaviour of the fit 
was investigated as a function of each of the input param¬ 
eters except seeing. Each parameter was varied over the 
range mentioned above, while keeping the other parame¬ 
ters fixed at the fiducial values. The results are shown in 
Figures and m For each parameter investigated, 12 

values over the range were chosen, and 40 realisations of 
each value were fitted for. 

The accuracy of the fit was relatively stable for Bind ^ 30° 
(Figure |^, and is stable for any input shear value in the 



8 (degrees) 

Figure 3. The error in the recovered shear values with varying 
inclination angle, • The black dashed line shows the fiducial 
inclination angle value used in the other fits. Symbols and input 
parameters as in Figure 1. 



Figure 4. The error in the recovered shear values with varying 
galaxy radius, p. The black dashed line shows the fiducial in¬ 
clination angle value used in the other fits. Symbols and input 
parameters as in Figure 1. 


range (—0.3 < 7 < 0.3), so long as the shear field is linear. 
The most notable change in sensitivity was with varying 
velocity signal-to-noise for Vmax/SV < 50 and with the 
number of pixels across the galaxy for p < 8 pixels. The 
results of this can be seen in Figures]^ and 

In each figure, the red dots show each individual realisation. 
The outer blue box shows the standard deviation, while the 
inner bars in the box show the uncertainty in the mean. 
Since the number of realisations per bin is reasonably large, 
the uncertainty on the mean is very small, such that it 
is hard to see in most bins in most figures. This suggests 
that we have performed enough realisations for the only 
uncertainties remaining in this suite of tests to be numerical 
ones. At worst these uncertainties contribute an error of 
~0(1%), which is an acceptable value. The exception to 
the above observation is the systematic bias introduced by 
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Figure 5. The error in the recovered shear values with varying 
noise, Yfrac- The black dashed line shows the fiducial inclination 
angle value used in the other fits. Symbols and input parameters 
as in Figure 1. 



Itrue (dimensionless) 


Figure 6. The error in the recovered shear values with vary¬ 
ing input shear, with seeing of s = 1, incorporated. Seeing has 
the effect of decreasing the magnitude of the fitted shear value, 
compared to the true shear. Symbols and input parameters as in 
Figure 1. 


seeing, which is discussed below. 

Once the effects of noise, shear strength, inclination angle 
and sampling resolution had been investigated, seeing was 
introduced, and its effect on the recovered shear accuracy 
was investigated by varying the input shear while holding 
the seeing fixed at s = 1 (corresponding to a spatial degra¬ 
dation kernel FWHM of 1 pixel), and then by varying the 
seeing while holding the input shear fixed at 7 = 0. The 
results are shown in Figures and [^respectively. 

The accuracy of the fit varies little with the size of the degra¬ 
dation kernel (Figure [^ for a fixed input shear of 7 = 0, 
however if the input shear is varied with fixed seeing, it be¬ 
comes apparent that seeing introduces a bias into the shear 
signal for any nonzero shear (Figure]^. This is because the 
seeing ‘smears’ the velocity map, washing out the shear sig- 



Figure 7. The error in the recovered shear values with varying 
seeing, s. Symbols input parameters as in Figure 1. 



Figure 8. The error in the recovered shear values with varying 
seeing (s), for four different input shear values in the model. Red 
corresponds to an input shear of 'yin = 0.3, magenta corresponds 
to an input shear of 7^^^ = 0.15, green corresponds to an input 
shear of 'jin = —0.15 and blue corresponds to an input shear of 
= —0.3. The boxes show the uncertainty of the mean (the 
inner horizontal lines), and standard deviation (the ends of the 
boxes), in each bin, while the circles show each individual realisa¬ 
tion. For each input value of 7, there are 40 realisations per bin, 
with 12 bins over a range s = (1,10). Fiducial values of p = 7.5 
pixels, Oinci = 35°, and Vmax/^V = 100 are used, as in the pre¬ 
vious models. The slight asymmetry in the plot about the x-axis 
is a consequence of the position angle (which is held constant in 
all of the realisations) of the model objects being neither 0° or 
±90°. The error in the fitted shear increases with shear size (and 
conversely is very small for small shears), so that a typical shear 
signal (7 < 0.1), when measured with good seeing, will have a 
very small degradation in the fit due to this effect. For example, 
for a true shear of 0.1 with a seeing of s = 1 pixel, the mean fitted 
shear is ~ 0.099, i.e. an error of 0.001. 
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Table 1. Guideline values for the investigated parameters, within which the DSM algorithm performs accurately. The “Range investi¬ 
gated” column shows the range in which each parameter was investigated. The “Lower limit” and “Upper limit” columns give the limits 
in which a realistic fit could be obtained, given the range investigated. The “Best fit” column gives the value of each parameter which 
resulted in the most accurate fit in the range investigated. 


Parameter 

Range investigated 

Lower limit 

Upper limit 

Best fit 

7 (unitless) 

-0.3 < 7 < 0.3 

-0.3 

0.3 

0.0 

dincl (degrees) 

10 < < 80 

30 

80 

45 

VmaxDV (unitless) 

1 < Vmaxl&V < 1000 

50 

none 

1000 

p (pixels) 

2 < p < 15 

7.5 

none 

15 

s (pixels) 

1 < s < 10 

none 

4 

1 


nal, and resulting in the measured shear having a smaller 
magnitude. Kaiser et al. 119951 have developed techniques 
to reduce the impact of such smearing in the case of 2D im¬ 
ages in traditional weak lensing applications by attempting 
to deconvolve the PSF from the image. It may be possible 
to implement similar techniques for velocity maps and 3D 
data in general, however we leave an investigation of this 
possibility to future work. 


To further investigate the effect of seeing on nonzero shears, 
a second set of model realisations were fit, over a range 
of degradation kernel sizes for four non-zero shear values; 
7 = -0.3,-0.15,0.15, and 0.3. As can be seen in Figure]^ 
the effect of seeing is to reduce the magnitude of the shear 
present in the data, and so increase the error. For a true 
shear which is positive (negative), increasing seeing give an 
increasingly negative (positive) error. This effect seems to 
saturate for very high values of s, corresponding to very 
poor seeing, and is relatively small for all realistic choices of 
seeing. For a seeing corresponding to a degradation kernel 
of s = 2.5, the error in the fit due to seeing is ~0(2.5%) 
at worst (i.e. in the case of the strongest possible shear). 
The error in the fitted shear increases with shear size (and 
conversely is very small for small shears), so that a typical 
shear signal (7 < 0 . 1 ), when measured with good seeing, 
will have a very small degradation in the fit due to this 
effect (< 0.1%). In light of this, this effect was not seen as 
being of immediate concern, so long as the seeing conditions 
for observing are reasonable. The slight positive bias in Fig¬ 
ure is due to the non-alignment of the projected distance 
between the lens and source and the position angle of the 
source, and is an expected effect. The uncertainty of the 
mean in this set of fits is larger than in the previous suite of 
tests because there are fewer realisations per bin in this case. 


Simulations were also undertaken to determine whether 
ellipticity in the point spread function (PSF) might bias 
the measurement of the shear using the DSM algorithm. 
A mock dataset was created, in which each object was 
convolved with a PSF with an ellipticity of e = 0.5. This is 
an extreme choice of PSF ellipticity, far greater than would 
be expected for any current telescopes. In spite of this 
however, the introduction of ellipticity to the PSF made 
a negligible (< 1 %) difference to the accuracy with which 
the DSM algorithm recovered the shears present in the 
mock dataset objects. Furthermore, no dependence between 
the PSF ellipse orientation and galaxy position angle was 
observed. 


Based on the errors per bin in the set of simulations de¬ 
scribed above, guideline acceptable fit ranges were specified 
for each of the input parameters. These define the parameter 
space in which the DSM method as currently implemented 
will be useful and effective. The guideline values for each 
parameter are presented in Table 


5 APPLICATION TO REAL DATA 

The requirements for the application of the DSM algorithm 
to a background galaxy are the following: 

• A reasonable expectation that weak lensing may be ap¬ 
plicable, vis an angular distance between the foreground 
and background galaxies that suggests an estimated shear 
> 0 . 01 . 

• The background galaxy has a suitable projection to 
the line-of-sight. Measurement of the shear from the veloc¬ 
ity map is not possible at inclination angles at or close to 
^inci = 0° and Sind = 90°. The measurement accuracy of the 
velocity map is minimum at these angles, and maximum at 
9 ~ 35°. This effect is taken into account in our investigation 
of velocity signal-to-noise in the previous section. 

• Ideally, the background galaxy is relatively isolated so 
that the disk motions are circular and undisturbed. 

• The projected separation between the lens (foreground) 
and source (background) galaxies does not coincide with ei¬ 
ther the major or minor axes of the source galaxy. 

• A measurable and regular rotation. In practice, this 
means that at high redshift there will be fewer suitable 
galaxies to which the DSM method can be applied. 

Non-circular motions in the disk of the background galaxy, 
such as tidal distortions or kinematic motions due to a bar, 
will result in well-defined signatures in the DSM fitting met¬ 
ric, and cannot be confused with the shear signal. In such 
cases it is not possible to fit for a small shear. 

In order to have confidence that DSM is fitting for true shear 
signals in lensed data, one needs to be reassured that where 
there is no lensing signal present, DSM recovers a null result. 
In the previous section DSM’s ability to accurately recover 
a null result in model data was demonstrated, however there 
are characteristics of real objects that were not included in 
the models. Examples of this include warp instabilities in 
the disks of galaxies, barred galaxies, ‘dumpiness’ and non- 
uniform structure within galactic disks, and turbulence in 
the outer regions of the galaxy (of particular relevance for 
radio data). 
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To assess how well a null result is recovered in un-lensed 
observational data, fits were attempted for a shear signal in 
the velocity maps of very low redshift galaxies. Since the 
probability of an object being lensed increases with distance 
(as the projected density of intervening matter increases), 
then the chance of a very low redshift galaxy being lensed 
is vanishingly small. With this in mind the shear was fit 
in a number of galaxies in the nearby universe, using data 
at radio wavelengths. Radio data was chosen for this first 
test since the models investigated in section [^have proper¬ 
ties analogous to that found in HI data cubes (specifically 
that the error maps are constant, rather than varying across 
the image). An investigation of optical integral field spectro¬ 
graph data is left for future work. 

Both ‘regular’ and ‘complex’ objects were fitted to. ‘Regu¬ 
lar’ objects are those that have a smooth velocity map with 
no warps, bars or clumps present. ‘Complex’ objects are 
those with warps apparent, barred galaxies, and those with 
clumpy structure and generally complex velocity profiles. 



RiglitAscension 


Data from The HI Nearby Galaxy Survey (THINGS) ( [mi¬ 
ter et al ?]|2008 1 was used, which provides very high angular 


and velocity resolution (7” and 5kms“^ respectively) 
observations of 34 objects with distances in the range 
3 <D< ISMpc (0.0007 < a < 0.00355). Since the velocity 
maps from THINGS are obtained from measurements in 
HI, they trace the HI abundance in the galaxies. The HI 
component of most galaxies extends out to much larger 
radii than the optical component, into regions where there 
is more turbulence and the disk is more likely to be warped. 
For this reason it is advantageous to mask this data beyond 
the extent of the optical component. Two objects were 
fitted from THINGS; NGC 3621 with morphological type 
SA(s)d, and NGC 5236 (M83) with morphological type 
SAB(s)c. The morphology of these two choices was of 
particular relevance, as it allowed for the effect of bars on 
the accuracy of the fit to be investigated. The HI emission 
and velocity map of NGC 3621 are shown in Figures [9 
and |10[ and those of NGC5236 are shown in Figures m 
and |12| In both cases the mask radius is 12 pixels and is 
shown by the black line. 


Two velocity maps were fitted from THINGS; one barred 
(NGC 5236) and one not barred (NGC 3621). Each were 
fitted with the same input parameters for the Monte-Carlo 
routine; 60 walkers, 200 burn in steps, followed by 300 
fitting steps. The fits for each object are presented in 
Tabled 


Initial guesses for the fit were obtained as follows: For the 
coordinates of the galaxy centre, the mean in the x and 
y positions of all nonzero pixels was used. For the central 
velocity, the mean of all nonzero pixels was taken. The 
position angle was estimated by eye. For all the fits an 
initial shear guess of 7 = 0 was used. So long as it was given 
sufficient time to converge, the algorithm was reasonably 
robust to poor initial guesses. With poor guesses up to 600 
burn in steps were required to obtain a fit of equivalent 
accuracy. Circular masks were added at approximately the 
extent of the regular velocity field, in the inner regions 
of the objects. This is roughly 1/3 of the radio extent for 
both NGC 3621 and NGC 5236. The use of circular masks. 


Figure 9. The zeroth moment (i.e the HI emission image) of 
NGC 3621 from THINGS. The mask radius is shown by the black 
line. The Fitted parametes for this galaxy are given in Tables 
and |3]for the masked and unmasked fits respectively. 



Figure 10. The first moment (i.e. the velocity map) of NGC 
3621 from THINGS, with the mask radius shown by the black 
line. The rotational velocities are based on the HI emission line. 
The Fitted parametes for this galaxy are given in Tables |^and|^ 
for the masked and unmasked fits respectively. 


necessitated to eliminate bias from complex extended 
structure, discourages fits involving strong distortions, and 
therefore damps the fitted shear values. This is an effect 
that needs to be addressed, however it is beyond the scope 
of the current analysis. 

As can be seen the accuracy with which the fits have recov¬ 
ered the expected null result is most affected by the appli¬ 
cation of a mask to the object. The warps and turbulence 
prevalent in the outer disks of even un-barred THINGS ob- 
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Table 2. Fitted parameters for the (masked) THINGS. Fits were obtained with a Monte-Carlo paramter space search using 60 walkers, 
200 burn-in steps and 300 iterations. 


Object 

7 

(unitless) 

^7 

(unitless) 

0 

(degrees) 

(degrees) 

3^0 

(pixels) 

<53^0 

(pixels) 

yo 

(pixels) 

(pixels) 

Vb 

(kms“^) 

5Vb 

(kms“^) 

NGC 3621 
NGC 5236 

-7.91 X 10“® 
5.22 X 10“® 

7.75 X fO-"^ 
1.85 X 10“® 

-9.65 

-119.70 

7.28 X 10"® 
1.46 X 10-1 

35.00 

41.03 

9.1 X 10“® 

1.2 X 10-® 

36.00 

40.97 

8.7 X 10“® 
1.2 X 10“® 

727.78 

509.56 

1.262 

1.298 
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Figure 11. The zeroth moment (i.e the HI emission image) of 
NGC 5236 from THINGS. The mask radius is shown by the black 
line. The Fitted parametes for this galaxy are given in Tables 
and for the masked and unmasked fits respectively. 
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Figure 12. The first moment (i.e. the velocity map) of NGC 
5236 from THINGS, with the mask radius shown by the black 
line. The Fitted parametes for this galaxy are given in Tables 
and [^for the masked and unmasked fits respectively. 



Table 3. Fitted parameters for the un-masked THINGS galaxies. 
Fits were obtained in the same manner as those presented in 
Table except that the objects were not masked. 


Parameter 

NGC 3621 

NGC 5236 

7 di ^7 (unitless) 

0 di (degrees) 

3^0 i Sxq (pixels) 
yo ± Syo (pixels) 

Vb i ^Vb (kms“^) 

0.049 ± 0.0057 
-9.66 ±0.415 
34.83 ± 0.069 
35.31 ± 0.021 
730.96 ± 0.580 

-0.234 ±0.2587 
-135.11 ± 7.414 
53.13 ± 13.217 

23.86 ± 13.513 
377.14 ± 91.497 


jects result in biasing of the fits, however masking to the 
optical radius of the object rectifies this. To illustrate the 
degree of improvement provided by this masking process, 
the THINGS objects were fitted for again without masking. 
The results of this are presented in Table It is clear from 
this table that masking to remove spurious signal from com¬ 
plex velocity structures, and retaining only the inner regular 
region of the velocity map is critical for obtaining good fits. 


6 CONCLUSIONS 

A new method for measuring gravitational shear from weak 
lensing has been described and tested. The method is based 
on the premise that symmetries in the velocity map are 
lost when it is lensed. DSM describes a process by which 
the deviation from symmetry might be measured and a lin¬ 
earised shear field determined. Systematic tests were per¬ 
formed to investigate the robustness of the algorithm, and 
requirements on the data quality for a high fidelity mea¬ 
surement. These tests also demonstrated the ability of the 
algorithm to recover a null result in synthetic data. The al¬ 
gorithm was then tested on images of nearby galaxies from 
the THINGS radio survey. A null result was obtained af¬ 
ter appropriate masking of the image. The results of these 
investigations, presented in Tables and have demon¬ 
strated that the DSM algorithm will recover a null result, 
both in simulated and real data. The range of parameters 
within which the DSM algorithm will be useful has also been 
established. These results are strong motivation to develop 
the DSM algorithm further. The first step has been made 
towards quantifying the errors and uncertainties expected 
when DSM is used to recover a shear signal from observa¬ 
tional data. Future work will investigate velocity maps from 
higher redshift galaxies that have a high probability of be¬ 
ing lensed to determine whether a non-zero shear signal is 
observed. This technique is most effective where the back¬ 
ground source is an isolated, undisturbed, stably-rotating 
disk. We have found that such systems not only exist but 
are reasonably easy to find (de Burgh-Day et al, in prep.) 
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in large datasets such as the Galaxy and Mass Assembly 
(GAMA) Survey ( Driver et al.|20rT Liske et al. in prep). 
An initial application will involve using this method to make 
measurements of individual dark matter halo masses, and 
measurements of cosmological parameters, including erg, and 
Ho. 
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